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We compute the frequency and magnetic field dependencies of the reflectivity R{uj) in layered super- 
conductors with two alternating intrinsic Josephson junctions with different critical current densities 
and quasiparticle conductivities for the electric field polarized along the c-axis. The parameter a de- 
scribing the electronic compressibility of the layers and the charge coupling of neighboring junctions 
was extracted for the SmLai-^Sr^CuCU-s superconductor from two independent optical measure- 
ments, the fit of the loss function L(oj) at zero magnetic field and the magnetic field dependence of 
the peak positions in L(u>). The experiments are consistent with a free electron value for a near the 
Josephson plasma frequencies. 



The Josephson plasma resonance (JPR) observed in 
the microwave absorption Q, and in optical reflection 
[§-|5| and transmission measurements HQ], has proven 
to be a powerful method to study the properties of 
highly anisotropic layered superconductors, such as vor- 
tex phases Ij]]^]. In particular, the spatial dispersion of 
the JPR, w p (k, k z ), parallel (k) and perpendicular {k z ) of 
the superconducting layers reflects the inductive coupling 
of junctions due to intralayer currents and the charge cou- 
pling due to the variations of the electrochemical poten- 
tial on the layers, respectively. Determining the latter is 
essential for understanding the coupled dynamics of the 
stack of intrinsic Josephson junctions in cuprate super- 
conductors, e.g. with respect to coherence in THz emis- 
sion p|JlO[| . From a fundamental point of view, it contains 
unique information about the electronic structure of the 
superconducting Cu02-layers, namely their compressibil- 
ity, which is hard to obtain otherwise. Beyond this, from 
the damping of the JPR the c-axis quasiparticle (QP) 
conductivity at the JPR frequency can be extracted. 

The parameter a characterizing the c-axis dispersion 
of the JPR in crystals with identical junctions is diffi- 
cult to observe in bulk microwave, transport fl(i| or opti- 
cal experiments, because mainly modes with small k z are 
excited, although for grazing incidence the JPR peak am- 
plitude can depend on a JHJ. The latter is irrelevant for 
the experiments in , which are performed with inci- 
dence parallel to the layers, which will be studied here. In 
the following we will show, how a can be extracted unam- 
biguously for layered superconductors with a superstruc- 
ture in c-direction from two independent measurements 
of the loss function and the magnetic field dependence of 
the plasma frequencies. 

Recently, these experiments have been performed on 
SmLai_ a: Sr a ;Cu04_5, where magnetic SHI2O2 and non- 
magnetic La2-xSr a; 02-5 layers alternate in the barri- 
ers between the CuC>2-layers ||-^]. In Refs. the 



reflectivity R and transmission of the electromagnetic 
wave propagating along the layers and with the elec- 
tric field along the c-axis (cf. Fig. la) were used 
to extract the effective dielectric function e e ff(u>) with 
the help of the Fresnel formulas, e.g. R{w) = [1 — 
\J e c ff (w)]/[l + yj e c ff (lo)]. The authors reported two 
peaks with quite similar widths in the loss function, 
L(u>) = Im [— l/e c ff(w)], and a quite large ratio of 
the peak intensities L(lo\) / L{u>2) between 10 and 20, see 
Fig. 2. In a two-junction model with two different plasma 
frequencies, but identical quasiparticle conductivities and 
no charge coupling the peak amplitudes in the loss func- 
tion are quite similar, see Fig. 2 at a — 0. This disagree- 
ment with the experimental data cannot be explained by 
a strong frequency dependence of the conductivities, as 
the widths of the two peaks are too similar. Instead, 
it was argued in |5|fT^] that the c-axis coupling plays a 
crucial role for the peak intensities. 
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FIG. 1. Geometry for measuring (a) the reflectivity with 
parallel incidence and (b) the microwave absorption in a cav- 
ity with the electric field polarized along the z-axis. 

However, in the derivation of e e ff in Ref. jFJj dissipa- 
tion was not introduced in the Maxwell equations, but 
arbitrarily in the final expression for e e g calculated with- 
out dissipation. In this Letter we correct these results 
accounting for different tunneling conductivities in the 
junctions in accordance with their different critical cur- 
rent densities. We compare the theoretical L(oj) with 
the experimental data in ||Jl| and extract at the JPR 
peaks the parameter a « 0.4, which is the free electron 
value. Further, we show that the microwave absorption 
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in the spatially uniform AC electric field applied along 
the c-axis (see Fig. lb) is also determined by L(ui). 

Recently, Pimenov et al. measured the dependence 
of the plasma frequencies on the magnetic field B || c. 
We show that the field dependence of the peak positions 
in L(u>) alone allows us to extract a 0.4, which is an 
independent confirmation of the fit of L(ui) for B = 0. 

We also show that the QP conductivities extracted at 
two different frequencies do not show the ^-dependence 
as anticipated for gapless d-wave pairing. We argue that 
the properties of the Cooper pair tunneling via the mag- 
netic Sm ions may overshadow the w-dependence of the 
QP conductivity. 
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FIG. 2. Experimental data [13] together with the calcu- 
lated loss function L(ui) (solid: a — 0.44, o\ = ra^ = 0.12, 
r — 0.38, ojo.i/c = 7.2 cm -1 , eo = 18; dot-dashed: a = 0, 
5-1 = r<? 2 = 0.12, r = 0.24, w ,i/c = 8.3 cm -1 , e = 13; 
dashed: a = 0, 02 = 15<7i = 0.2, r = 0.24, Wo,i/c = 8 cm -1 , 
eo = 15). 

We consider the crystal with two alternating Joseph- 
son junctions characterized by different critical current 
densities, Ji, and the c-axis tunneling conductivities, 07, 
due to different tunneling matrix elements. We assume 
that all other parameters of the junctions are identi- 
cal, as their distinction would lead to negligible correc- 
tions in the following. Without the charging effect the 
c-axis bare plasma frequencies Wo,i are related to J; as 
u>q ; = 8ir 2 cs Ji/ 6q$ 0, where $0 is the flux quantum, eo is 
the high frequency c-axis dielectric constant and s is the 
interlayer distance. We neglect nonequilibrium effects 
p0| , p^t in the distribution function of the quasiparticles 
assuming that the plasma frequencies are well below the 
charge imbalance and energy relaxation rates. 

To find the reflectivity R(uj) in parallel incidence and 
the microwave absorption, we use the Maxwell equations 
inside the crystals accounting for supercurrents inside the 
2D layers at z = ms and interlayer Josephson and quasi- 
particle currents determined by the difference of the elec- 
trochemical potentials in neighboring layers: 
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Here p m is the chemical potential in the layer m, V mjTn +i 
is the difference of the electrochemical potentials, Lu a0 — 
c/X a b^ye a o is the in-plane plasma frequency and e a o is the 
high frequency in-plane dielectric constant. The function 
/ is defined as f m ,m+i(z) = 1 at ms < z < (m + l)s and 
zero outside this interval. To obtain Eq. (|^) for small am- 
plitude oscillations we expressed the supercurrent density 

(s) 

J m,m+i = Jl sin (p m>m+ i ~ Ji<p m ,m+i via the phase dif- 
ference ifim,m+i — 2iV mjm +i/Hu). Further, Co 2 = uj 2 (1 — 
iAiraiU) / u)q ; co) _1 takes into account the dissipation due 

to QP tunneling currents, J^fm+i = <T iKn,m+i/ es j which 
are depend on the different conductivities 07 in the junc- 
tions I — 1,2 and the difference V m . m +\ of the electro- 
chemical potentials. We express the difference of the 
chemical potentials p m via the difference of the 2D charge 
densities, p m , as p m - p m+ i = (4irsa/e )(p m ~ p m +i), 
where the parameter a — {cq/ Aires) {dp/ dp) character- 
izes the interlayer coupling. In the model of 2D free elec- 
trons we find dp/ dp = irH 2 /(em*). For an effective mass 



to* w 1 — 2m e , as expected from ARPES [jlf 
and eo = 20 we estimate the order of a as « 
The solution of the Eqs. (0)-(|) is 
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c m exp(igz) - dm exp(-igz)}, 
d m exp(-igz)] + aA m , (5) 



where now < z < s in each junction, g 2 = uj 2 e a o/(c 2 a) 
and a -1 = 1 — c 2 k 2 /eou> 2 . The Maxwell boundary condi- 
tions for By and E x across the layers lead to a set of equa- 
tions for c m , d m and A m , which we solve using the small 
parameters b = gs/2 ~ s/X c <C 1 and (3 — s 2 /2X 2 b a <C 1. 
In order to find the reflectivity for parallel incidence, cf. 
Fig. la, we consider the case sk z <C b < (3, which is 
fulfilled for incident angles S « 1, and obtain the dis- 
persion relation k x (uj) for waves propagating inside the 
crystal (for details see fill): 



e e ff(w) 



r(w — vi)(w — V2) + iS 
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rw 2 - (1 + r)(2a + l/2)w + iSi 



77^- (6) 



vi,2 = (1 + r)(l + 2a)(l T ^/T^p)/2r, 

4r(l + 4a) 
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Si = w^ 2 r(2a + l/2){a 1 + d 2 ), (9) 
5 = w 1/2 [{2a + l)rw(ai + a 2 ) - (1 + 4a)(cri + a 2 r)], 

where vi = u/f/uj^, u>i are the JPR frequencies, w = 
uj 2 /cJq V ai = 47Tf7;/e wo,i and r = uj 2 a /uj 2 2 < 1. 
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FIG. 3. Ratio of the peak amplitudes L(ui\) / X(wa) in the 
loss function depending on the squared ratio aii/wf of the 
JPR frequencies for different a — 0, ... , 0.65 (top to bottom 
plot) at <Ti = ra 2 . 
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FIG. 4. Ratio r = Wo,i/ w o,2 of the squared bare plasma 
frequencies depending of the squared ratio io\ /w| of the JPR 
resonances in L(uj) for different a = 0, .., 0.6 (from below). 

The reflectivity coefficient is given by the usual Fresnel 
expression R = (1 — ck x /u)/(l + ck x /u). Note that for 
nonzero [3 we obtain in principle two propagating modes 
inside the crystal, but the second mode gives a negligible 
contribution near the JPR peaks at j3 <ti 1, see Jl^] . The 
loss function is 

L{w) = Im [-l/e cff H] = (w 3/2 /2e ) x (10) 
a±(wr — 4a — I) 2 + <T 2 r 2 (u) - 4a - l) 2 
[r(w — vi)(w — v 2 )] 2 + S 2 
and shows resonances at the two transverse plasma 



bands, while the peak in Im[e e jf(o;)] is at vt — (1 + 
l/r)(2a+l/2). 

In s-wave superconductors with nonmagnetic ions in 
the barrier between the layers both the QP conduc- 
tivities 07 and the critical current densities J/ vary 
with I only due to their proportionality to the squared 
matrix element for the interlayer tunneling in agree- 
ment with the doping dependence <j c (x) oc w 2 (x) found 
for La2- a; Sr :c Cu04 p7| . We take the same relation, 
o\(yj\) I oi(u>2) = r for SmLai^Sr^CuC^-^, ignoring the 
frequency dependence of ai due to the d-w&ve pairing 
p8[ , which decreases cti(wi)/ct2(w2) below r. On the 
other hand, we also neglect that due to the magnetism 
of the Sm ions the critical current density J\ a j 
is suppressed, while oi is not affected, which suggests 
<7i/(72 > t without d-wave pairing ||l9| . The fact that our 
assumption o\ = ra-i turns out to be consistent with the 
experimental data suggests that both effects compensate 
each other. 

In Fig. H we present the dependence of the ratio 
of the peak amplitudes in the loss function L(uj), i.e. 
L(u)i)/ L(u>2), vs. the ratio o;f/w| for different values of 
a. Hence, this figure allows us to obtain the parameter 
a from the positions lui.2 of the JPR peaks and their ra- 
tio of amplitudes. Fig. || allows us to find the ratio of 
the squared bare frequencies, r = u/q 1/^02 f° r a ratio 
lui/uj2 when the parameter a is obtained from Fig. |^. 
The comparison of L(uj) in Fig. ^| with the experimental 
data from Ref. |L3| gives the best fit for a\ — r&2 = 0.12 
and eo = 18, a = 0.44, which is of the same order as the 
theoretical estimate for free electrons, cf. Tab. I for a fit 
of pP^,0. In the case a — a, high ratio o^/ai ~ 15 of 
the conductivities is necessary to reproduce the ratio of 
the amplitudes, which fails to describe the shape of the 
resonance at L02 correctly. 

The loss function L(lo) also determines the microwave 
absorption of a crystal in a capacitor, which induces a 
uniform (sk z — > 0, k x — 0) AC electric field 2£ cos(ojt) 
above and below the crystal, see Fig. lb. The crystal ex- 
citations are longitudinal in this case, i.e. E x = B y = 
in Eqs. (|l|)-(§). These equations together with the Pois- 
son equations near the top and bottom layers determine 
for large N the microwave absorption as 



V{oj) ^N- 1 <ri\V m , m +i/esf 
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Next we consider the dependence of the JPR peaks in 
L(u>) on the c-axis magnetic field B, which allows one to 
estimate the parameter a, without relying on the abso- 
lute amplitude of the spectra. For Josephson junctions 
B suppresses the critical current densities and the bare 
plasma frequencies, u>qi(B), and also broadens the res- 
onance peaks due to the formation of pancake vortices 
randomly misaligned along the c-axis [H. From Eq. (|?]) 
we see that the ratio of the resonance frequencies, 
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^jui = (i - x/i^va + v^p), 



(12) 



depends only on the quantity p given by Eq. (jq), 
which varies with i? only via the function r{B) = 
w o i(B)/^q 2 (B). When the functional form of ujq^{B) is 
known theoretically and p(B) is obtained from the mea- 
surement of W1/W2 via Eq. (|l2|), one can fit unknown 
parameters in loqj(B) and obtain a. It is pointed out 
that extracting a from L(oj, B) by fitting all parameters 
without further theoretical input is difficult due to the 
dependence of the effective c-axis conductivity on B || . 

The dependence loqj(B) in the decoupled vortex liquid 
H is known for crystals with equivalent junctions and can 
be used for alternating junctions as well, provided that 
the vortices in different layers are decoupled. This is sug- 
gested by huJo.i <S kT c in SmLao.sSro^CuO^a (cf. Tab. 
I) indicating Josephson coupling of the layers. Thus, for 
the vortex liquid phase with decoupled pancakes in high 
magnetic fields B > Bj = >$> Q \ 2 ab /{\ c s) 2 - 1 T the de- 
pendence of the plasma frequencies on B is M 



32Btt 3 c 2 sT' 
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This gives r(B 3> Bj) = r 2 (B = 0). Hence, we find 



P(B > Bj 



4r 2 (0)(l +4a) 
[l + r 2 (0)] 2 (l + 2a) 2 
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while for p(0) we use Eq. (g) with r = r(0). 

Consequently, the parameters a and r(0) can be found 
analytically without fitting from p(0) and p(B ^> Bj) ob- 
tained from JPR frequency measurements using the data 
of Ref. U . The dependence uj 2 ,oj 2 oc 1/B at fields above 
Bj ss 1 T shows that in this field range the decoupled 
pancake liquid is present. We use the reported field de- 
pendence of the peak frequency in Im[e ff(w)] for uj 2 (B), 
as the authors noted that they are quite close and the dif- 
ference will finally turn out to be ss 10 %. Then we obtain 
(for T < T c ) tul/uj 2 w 0.31 at B = and 0.21 at B > Bj 
or p(Q) ss 0.73 and p(B 3> Bj) ss 0.57 respectively. From 
Eqs. (||) and we obtain r(0) ps 0.55, a ps 0.4 and the 
bare frequencies are u>o,i/c — 6.6 cm" 1 and loq^/c = 8.9 
cm" 1 , cf. Tab. I. This estimate gives the parameters a 
and r(0) similar to those obtained above from the fit of 
L(oj) at B = and confirms independently the relevance 
of the c-axis coupling a and our assumption &\ ss rdi. 
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TABLE I. Extracted parameters for SmLai-^Sr^CuO^ 



In conclusion, we calculated the effective dielectric 
function e e ff, Eq. (||), for alternating junctions with 
charge coupling a and frequency independent QP con- 
ductivities <ti tx Lu 2 i . This allows to describe satisfac- 
tory the optical properties of the SmLao.sSro^CuOzi-s 
superconductor near the plasma frequencies. The pa- 
rameter a ~ 0.4 was extracted independently from the 
magnetic field dependence of the positions of the JPR 
peaks bJ\^{B) and from the fit of the loss function L(u>) 
at B = 0. Its value is expected to be universal in the 
cuprates and corresponds to the free electron value of 
the electronic compressibility, which can differ from the 
renormalized one particle density of states both for Fermi 
and non Fermi liquids |2(i]. We also show that the ex- 
tracted conductivities c; differ according to the differ- 
ent tunneling matrix elements, ai/a 2 ~ w o 1/^02- This 
relation is anticipated for tunneling between s-wave su- 
perconductors via nonmagnetic ions. Its fulfillment in 
SmLao.8Sro.2Cu04_5 might indicate that the gapless d- 
wave pairing in the layers reduces u\{uji) / 02^2) in a 
similar way as the magnetism of the Sm ions decreases 
Wg i/u}q 2 - This could be studied further by measuring the 
gap in the QP spectrum in the I-V curve as in Ref. [Q . 
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